Theory for a dissipative droplet soliton excited by a spin torque nanocontact 
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A novel type of solitary wave is predicted to form in spin torque oscillators when the free layer 
has a sufficiently large perpendicular anisotropy. In this structure, which is a dissipative version 
of the conservative droplet soliton originally studied in 1977 by Ivanov and Kosevich, spin torque 
counteracts the damping that would otherwise destroy the mode. Asymptotic methods are used to 
derive conditions on perpendicular anisotropy strength and applied current under which a dissipative 
droplet can be nucleated and sustained. Numerical methods are used to confirm the stability of the 
droplet against various perturbations that are likely in experiments, including tilting of the applied 
field, non-zero spin torque asymmetry, and non-trivial Oersted fields. Under certain conditions, the 
droplet experiences a drift instability in which it propagates away from the nanocontact and is then 
destroyed by damping. 

PACS numbers: 75.40.Gb 85.75.-d 75.40.Mg 76.50.+g 75.30.Ds 75.75.+a 
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I. INTRODUCTION 

The concept of a soliton as a localized particle-like 
excitation that preserves its shape can be extended to 
systems that are far from thermodynamic equilibrium 
through the concept of a "dissipative soliton"—. This 
allows us to analyze a broad range of physical, chemi- 
cal, and biological nonlinear systems in which localized 
excitations are observed. Driven magnetic systems, es- 
pecially those of technological interest, exhibit strongly 
nonlinear dynamics and are an ideal experimental do- 
main for exploring the dissipative soliton model. 

In this paper, we provide an analytical theory for a 
novel, localized oscillation mode in a spin torque oscil- 
lator with a free layer having perpendicular magnetic 
anisotropy. The salient features of this mode include 
a frequency well below that of uniform ferromagnetic 
resonance, a weak dependence of frequency on bias cur- 
rent, and a precession angle at the maximal value of 90°. 
Combining numerical micromagnctic simulations with an 
asymptotic analysis of the equations of motion, we iden- 
tify this mode as a dynamic, dissipative magnetic soliton 
that is closely related to the "magnon droplet" predicted 
by Ivanov and Kosevich in 1$772£. The mode central 
region exhibits magnetization pointing nearly opposite 
to its equilibrium direction with a perimeter manifesting 
90° precession. From our asymptotic analysis, we derive 
conditions on perpendicular anisotropy and bias current 
for the nuclcation and existence of the dissipative droplet. 
Using our numerical simulations, we analyze the stability 
of the dissipative droplet soliton as a function of applied 
magnetic field, bias current, and spin torque asymmetry. 

Solitons in conservative systems occur when nonlinear 
terms in the equation of motion balance the effects of 
dispersion^. A classic example is a light pulse moving 
in a lossless optical fiber: the change of refractive in- 
dex with frequency (dispersion) tends to make the pulse 
spread out, but for a certain pulse shape, the change 
of refractive index with light intensity due to the opti- 



cal Kerr effect (nonlincarity) exactly balances the dis- 
persion. Pulses having this shape can propagate without 
spreading and are called solitary waves or solitons. The 
balance between nonlincarity and dispersion typically al- 
lows for the existence of a continuous family of solitons 
that can be excited in the system, rather than a single 
solution. In the optical fiber example, the family can 
be parametrized, for example, by pulse amplitude, and 
there is a continuous range of amplitudes that satisfy the 
soliton balancing condition. 

Dissipative solitons^ are characterized by an additional 
balancing condition between gain and loss that typically 
allows only a single solution for a given set of external 
parameters. Although conservative soliton models can 
explain weakly nonlinear behavior seen in magnetic sys- 
tems of exceptionally low damping^, damping is not a 
small effect for many magnetic systems of both funda- 
mental and technological interest. By combining clas- 
sical soliton theory with bifurcation theory of nonlinear 
dynamics and concepts of self-organization^, the dissipa- 
tive soliton concept provides a framework for describing 
a broad range of soliton-likc behaviors. Here, we ap- 
ply this concept to a nanoscale ferromagnetic system in 
which both damping and a driving force (spin torque) 
are important. 

Spin torque^— occurs when a current is driven through 
a structure with alternating magnetic and nonmagnetic 
layers in which spin-dependent conductance at the in- 
terfaces results in a spin-polarized electron flow. When 
the polarized electrons enter a ferromagnetic layer whose 
magnetization M is not collincar with the electron spins, 
the transmitted spins are rotated toward M and the an- 
gular momentum absorbed by the ferromagnetic layer is 
known as the spin torque. Typical devices have two ferro- 
magnetic layers through which current is driven: a thick 
"fixed" layer that determines the direction of electron 
polarization, and a thin "free" layer whose orientation 
can be readily changed by spin torque. For current of 
the appropriate polarity, spin torque opposes the intrin- 
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sic damping torque in the system, and currents above 
a threshold produce dynamic states in which M of the 
free layer can be manipulated without applying a mag- 
netic field. This effect has been used to control switching 
of nanoscale magnetic elements^, with potential appli- 
cations in computer memory and data storage. The ef- 
fect has also been used to produce coherent, frequency- 
tunable microwave oscillations^ in a nanoscale device 
known as a spin torque oscillator (STO), with poten- 
tial applications in integrated microwave circuits for mix- 
ing and active phase control. Recent reviews cover the 
physics of spin torque^ and its possible applications^. 

The equations of motion for M in the presence of 
spin torque (presented below) are inherently nonlinear, 
and their full solution for a general case is often stud- 
ied by use of numerical methods. Analytical methods 
can sometimes be applied by invoking restrictions such 
as high symmetry, spatially uniform M (the "macrospin" 
model) , and small precession amplitude (small angle be- 
tween M and its equilibrium direction) . These restricted 
cases have been used to explain experimental results with 
mixed success. The local nature of spin torque allows it 
to drive large amplitude excitations in which M varies 
on the scale of the magnetic exchange length (typically a 
few nanometers), something that applied magnetic fields 
cannot do. As we show here, this regime of strongly non- 
linear, strongly nonuniform, sustained magnetodynamics 
is amenable to theoretical examination using numerical 
and analytical approaches. This regime is also experi- 
mentally accessible in STOs. 

We note that a different type of magnetic soliton gen- 
erated by spin torque, called a spin wave "bullet" , was 
predicted by Slavin and Tibcrkcvich in 2005 to occur 
in the point-contact geometry with magnetic films ex- 
hibiting in-plane oriented anisotropy and in-plane applied 
magnetic fielcU^. For this case, the precession frequency 
decreases with increasing current, which can result in lo- 
calization if the frequency falls below the bottom of the 
spin wave band at the ferromagnetic resonance (FMR) 
frequency. The weakly nonlinear bullet soliton is a solu- 
tion to a Nonlinear Schrodinger type equation with third- 
order nonlinearity in the excitation amplitude. As such, 
its predicted experimental signature consists of subtle 
shifts in microwave output frequency and threshold cur- 
rent relative to that expected for a non-localized mode. 
In contrast, the droplet soliton studied here exhibits dra- 
matic differences in behavior from that of a non-localized 
mode. This is due to the fact that it is a strongly nonlin- 
ear solution of the full equations of motion, rather than 
simply a third order expansion. 

Domain wallsi^, magnetic bubbles^, and vortices^! 
are examples of well studied, strongly nonlinear, localized 
structures that occur in magnetic materials. The droplet 
differs from these static structures in that it is inher- 
ently dynamic; the frequency of spin precession within 
the droplet is always greater than zero. In this work, 
we focus on the two-dimensional (2D), non-topological 
droplet, but we note that droplets in two and three- 



dimensions come in topological flavors as well^. 

We begin in the next section by presenting an asymp- 
totic analysis of the model equations for the dissipative 
droplet in a high-symmetry case. We will also derive 
the droplet's frequency vs. current relation in this sec- 
tion. Section IIIII is devoted to the study of droplets in 
physically realistic situations incorporating the current- 
induced Oersted field as well as canting of the applied 
field and fixed layer. Section IIVI details experimentally 
accessible nucleation conditions for a droplet that take 
advantage of a small amplitude instability. In Sec. [Vj we 
discuss possible extensions of the theory and we relate 
the droplet to other excitations in thin magnetic films. 
We conclude in Sec. I VII with a summary of the droplet's 
unique properties. Appendices |A1 andlBlprovidc details of 
our stability calculation and numerical method, respec- 
tively. 

II. DROPLET IN A NANOCONTACT 

We consider the Landau-Lifshitz-Slonczewski equation 
in non-dimensional formic describing the free layer 
magnetization in polar coordinates (m = m(p, ip, t); bold 
symbols represent vectors in M 3 or R 2 , contextually de- 
pendent; e.g., m = [m x ,my,m z }) 

— = m x heff-am x (m x h cff ) 

precession damping 

+ aV{p) m x (m x m f ), m : R 2 — > § 2 , (1) 
1 + vra • mf 

V v ' 

spin torque 

h c ff = V 2 m + h + h oc + (hk - l)m z z, |m| = 1. 

Figure [1] is a schematic of the nanocontact geometry and 
coordinate systems considered in this work. The magne- 
tization m = M/A/g and fields (e.g., ho = Hq/M s ) are 
normalized by the saturation magnetization M s , time is 
normalized by 2tt /\j\fioM s (7 is the gyromagnetic ratio 
and /io is the permeability of free space), and space is nor- 
malized by the exchange length ^ cx = */ D/\j\noM s h (D 
is the exchange parameter and h is the reduced Planck's 
constant). We consider a free layer whose thickness 
5 is much smaller than the magnetic excitation wave- 
lengths, so that a 2D model with local dipolar fields is 
justified^. The precessional term is driven by the effec- 
tive field h e ff incorporating exchange V 2 m, a uniform ap- 
plied field ho = /io[sin6*o, 0, cos^o] (ho > 0), the current- 
induced Oersted field h oc = —f oc (p)tp (the definition of 
f oc is given in Ref. [l^ and summarized in Sec. IIII Al) . 
the perpendicular field hkm z z due to crystal anisotropy 
(/ik > 1), and the demagnetizing field for a thin film 
— m z z. The Landau-Lifshitz damping coefficient is a. 
The spin torque term involves a = I/I c , where / is the 
current, I c = (A 2 + l)2M 2 e^ Trr 2 J/(hPX 2 ) [P and A > 1 
are the polarization and spin torque asymmetry param- 
eters, respectively^, e is the electron charge, is the 
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In the symmetric case of Eq. ([2]), we can remove the 
explicit dependence on h k from Eqs. and (UJ) by in- 
troducing the following rescaling: 



(p 











7\tn 














J ^ 



p = p'/y/h k -l, t = t'/{h k -l), 

a = (hk - l)a', h = (h k - l)h' Q , p* = p'J y/h k -l. 

(5) 

Recall that we are assuming h\n > 1. Then, with the 
scalings in Eq. ([S]) and dividing Eqs. and (@| by /ik— 1, 
we can, without loss of generality, take — 1 = 1. For 
the rest of this Section [Til we will use the scalings in 
Eq. ([5]) so that hu. -> 2. 

In this Section, we consider localized magnetization 
configurations that satisfy 



FIG. 1: Schematic of nanocontact (yellow disk of radius p*), 
the domain of m (R 2 in polar coordinates with radius p, az- 
imuthal angle ip), the range of m (unit sphere 8 2 with polar 
angle O, azimuthal angle <£>), orientations of the applied field 
ho and fixed layer mf. 



lim Q(p',<p,t') = 0, 
p' — >ao 

with sufficiently rapid decay. As such, we define the mag- 
netic energy in terms of exchange and anisotropy energy 
via 



nanocontact radius), v — (A 2 — 1)/(A 2 + 1) (0 < v < 1), 
V(p) = H(p* — p), where H is the Heaviside step func- 
tion, p* = r*/^ G x is the nondimensional contact radius, 
and rrif = [sin 8f, 0, cos 9f] is the uniform fixed layer mag- 
netization. 

Numerical computations of Eq. ([TJ will be presented 
in Sec. IHII bv means of a method discussed in Appendix 
|B1 For the rest of this section, we will focus on a high- 
symmetry case where 



high symmetry case: 8q = 0, h Q 



0. 



(2) 



For the analysis of this case, it is helpful to 

consider Eq. (TTJ) in spherical coordinates m = 
[cos <f> sin 0, sin $ sin 0, cos 0] (sec Fig. [T]) 

^ = F[Q, *] - aG[Q, $] + <jV(p)P e [Q, $], (3) 

sine— = G[Q,$] + aF[Q,$\ + aV(p)P*[Q,§], (4) 



where 

F[0,$] 
G[0,$] 



sin 0V 2 $ + 2 cos 0V$ • V0, 
1 



-V 2 



ho sin . 



sin20(|V$| 



1) 



— cos cos <f> sin Of + sin cos Or 

Pq 10 $1 = , 

1 -I- ^(cos <& sin sin Of + cos cos Of ) ' 

p r@ $1 _ sin $ sin 6 f 

* ' 1 + j/(cos<I>sm0sin#f + cos0cos6*f) 

The polar angle satisfies < < n, while the az- 
imuthal angle $ is interpreted modulo 2-7T. 



£[©,$] = i 



|V0| 2 + sin 2 0|V$| 2 + sm^0 

anisotropy 



exchange 



Note that the damping and spin torque terms break en- 
ergy conservation 



d£ 
~dt 



[6,*] 



/ 1 




/ffi 2 





(G[0, $] - sin0)P [0, $] - ^[6, *]P*[0, $] 



(G[0, $] - ^ sin 0) 2 + F[0, $] : 



dr'. 



(0) 



A. Conservative Droplet Soliton 



In the absence of damping and spin torque (a = 0, 
a 1 =0), Eqs. © and Q admit a continuous family of 
exponentially localized, non-topological solutions known 
as magnon droplet solitons 2 -^. These solutions can be 
parametrized by the frequency ujq and written as 



= o (p;o;o), 
$ = (w + h' a )t', 

= [wo(hk - 1) + ho\t. 



They satisfy a balance between exchange (dispersion) 
and anisotropy (nonlinearity) through the nonlinear 
eigenvalue problem F[0o, (w + h' )t'] = 0' G[0o, (wo + 
h' )t'} = (uio + ft-o)sin0o or 

(|^ + 7^) eo ^ sin2eo + CJoSineo ^ ' (8) 
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FIG. 2: Conservative droplet profiles with m z — cos0q. 



with the boundary conditions 



0, lim e (p';w ) =0. (9) 

p' — >oo 



The conservative droplet solutions in Eqs. ([7]), ©, and 
© have the following properties. The polar angle 0o 
varies with radial distance and is independent of time; 
thus the spatial distribution of m z = cos Oq is static 
and azimuthally symmetric. The azimuthal angle <£> is 
independent of position and varies linearly in time; thus 
all points precess at a common frequency and in phase. 
It has been shown thalj 2 -^ ujq satisfies 



< uj < 1. 



(10) 



Therefore, the total precessional frequency ujq + h' varies 
between the Zeeman frequency h' and the frequency of 
spatially uniform precession about = 0, 1 + h' Q (the 
FMR frequency). 

While the conservative droplet does not have a closed- 
form analytical expression, we calculate it by numeri- 
cally solving Eq. ([8]) subject to the boundary conditions 
© with the function bvp4c in Matlab®. A plot of sev- 
eral conservative droplet profiles is shown in Fig. [2] We 
see that the amplitude at the origin 1 — m z (0;ujQ) = 
1 — cos Oo(0; wo) decreases as the frequency is increased. 
While it may appear in Fig.[5]that m z (0;wo = 0.15) = 0, 
in fact all conservative droplets with < ujq < 1 satisfy 
m z {Q\Lo Q ) > —1 owing to their non-topological structure. 

The energy for the droplet £ = £{<d ,(uj Q + h' )t'] 
satisfies 2 -^ 



dij 



< 0. 



(11) 



The fact that the energy is a decreasing function of fre- 
quency has been used to argue that the conservative 
droplet is stable in a 2D material^. The 2D conservative 
droplet embedded in an infinite, 3D magnet is known to 
be unstable. However, preliminary studies suggest that, 
for sufficiently thin films, the 2D conservative droplet is 
stabilized. This work is beyond the scope of this paper 
and will appear elsewhere. 



As we now show, the conservative droplet soliton 
can be generalized to damped/driven systems such as a 
nanocontact. Whenever we refer to "droplet", we mean 
the dissipative droplet studied in the future sections. We 
will always use "conservative droplet" to describe the so- 
lution of Eq. (|8]) that is monotonically decaying to zero 
as p — > co . 



B. Dissipative Droplet Soliton 

We now extend the analysis of Kosevich, Ivanov, and 
Kovalcv^, to the case of the dissipative droplet solution, 
where damping is no longer assumed to be negligible, 
and spin torque is included in the analysis. In addition 
to the balance between exchange and anisotropy that was 
required in the conservative droplet case, a further bal- 
ance between uniform damping (loss) and spatially local- 
ized spin torque (gain) will be derived for the droplet to 
be sustained. We will assume that the spin torque and 
damping, while small, are not zero and are of the same 
magnitude so that 

(7 = 0(a), < a « 1. 

We look for an asymptotic solution in the following form 

6(r',<') = Q (p';u J )+aQ 1 (r',t') + --- , 



$(r',t') = (w + h' )t' + a- 



$i(r',Q 
sine (p'; uj) 



(12) 



We have set uio — > uj to distinguish the frequency of 
the droplet u> from that of the conservative droplet u)q. 
We will conclude that there is no frequency shift due to 
damping and spin torque so that ui = uiq here. However, 
other perturbations beyond those considered here could 
lead to a frequency shift. Inserting the expansions (fT2)) 
into Eqs. [31 21 and equating like orders in a gives the 
following equations for the perturbations Oi and $i: 



— = -(uj + ti )sm@ 



dtf 



(13) 



+ -V(p')Pe[e ,(uj + h' )t' 



0*i 



dt' 



a 
a 

a 



y + L e Oi = V(p')Pz[Q , {uj + h' )t% (14) 



where the sclf-adjoint, Schrodingcr operators L$ and Lq 
are 

L* = - S ^[Q ,(uj + h' )t']^— 
09 smBo 

= -V' 2 -4 +cos 2 9 - wcos9 , 

dp' 

SG 

Lq = ^q[ o, (uj + h' )t'] - (uj + /i )sin9 
= -V' 2 + cos 20 o - u;cos6 . 
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Note the following important property: 

L<£,sin0o = O. (15) 

The rest of this section is concerned with the solu- 
tion of the pcrturbative equations (|T5t and ((H)) in two 
separate cases. We use standard solvability arguments 
for forced, linear differential equations to remove secular 
terms (sec e.g. Ref. H3) in order to derive an expression 
for the current a' at which the balancing condition for 
the droplet is achieved. 



1. Case Of = 

We first consider the case where the fixed layer is ori- 
ented perfectly normal to the film plane so that 6f = 0. 
In this regime, we can study the effect of variable spin 
torque asymmetry v. 

In addition to the scalings in Eq. , we can also scale 
out the dependence on the applied field h' when Of = 
with the substitution 




-15 -10 -5 5 10 15 

x' 



FIG. 3: Dissipative droplet soliton for the high-symmetry case 
with Of =0. The color scale corresponds to m z , while the 
vector field corresponds to the in-plane components (m x ,m y ). 
The circle here and in future plots represents the boundary 
of the nanocontact. Parameters are cr sus /a = 0.94, u = 0.17, 
pi = 5.24, and v = 0. 



$ = $' + h' t'. 



(16) 



Then, 6 and satisfy Eqs. ^ and (gjl, as before, but 
with h' = 0. This transformation shows that the vertical 
applied field simply shifts the precessional frequency by 

K- 

We seek a solution to Eqs. (|T5|) and p^|) in the form 
f = 0, h' = 0, 0i = 0, and $i(r',i') = (i.e., 
linear phase evolution with time) which results in the 
following non- homogeneous equation for $^1 



i^^j = — wsinGo 
$[(0) = 



<x'VV)sin0 o 
a(l + ^cos0o) 



(17) 



This equation is solvable if and only if the nonhomoge- 
neous terms are orthogonal to the kernel of By use 
of Eq. (fl"5jl , multiplying the right-hand side of Eq. (fl~7|) 
by p sin O and integrating from to infinity we ob- 
tain the existence condition for a dissipative droplet, 
c' = cr sus (u;), where 



c S us(w) = auj 



r o °°sin 2 Q (p';u})p'dp' 
j V(p') S m 2 Q (p';uj) ■ 
J 1 + i> cos e {p';uj) H 



(18) 



The choice a' = <r sus (w) singles out a specific value for 
the current as a function of the droplet frequency < 
lu < 1. This value of the current provides the necessary 
balance between damping and spin torque, in addition to 
the balance between exchange and anisotropy, in order 
to sustain the droplet. Therefore, we call tr sus (w) the 
sustaining current. 



We can also understand the choice for the sustaining 
current (|18p by computing the rate of change in the en- 
ergy for the droplet from Eq. (J6j) : 



— [0 O) u;i'+a$i/sin0o] 
at 

f sus(L0)V(p') 



sin 2 0n 



aui sin 2 O 



1 + z/cos 0o 
dr' + 0{a 2 ) = 0{a 2 ). 



(19) 



Therefore, the energy of the droplet is approximately con- 
served for the choice a = a sus . Figure[3]shows a snapshot 
in time of a dissipative droplet. 

Figure 3] represents the numerical evaluation of 
Eq. (fT8|) and shows the dependence of er sus on uj. The 
droplet frequency has two branches as the sustaining 
current is varied, but only one branch is stable. For a 
given droplet frequency uj, consider increasing the current 
slightly above the sustaining value, a' = a sus (io) + 5a'. 
From Eq. (|19[) . the droplet energy will increase, and from 
Eq. (fPTj) , this increase in energy corresponds to a decrease 
in droplet frequency for a stable droplet^. Thus, the up- 
per branch in Fig. |4l for which increasing current causes 
an increase in frequency, is unstable. 

Figure 2] shows that there is a minimum sustaining 
current er™ 1 ", maximum frequency o; max , and minimum 
m z in the center of the droplet m™ ax where 

o™?= min o sns {uj) = a sus (^ max ), 

"6(0,1) 

m I 2 nax = cos0 o (O;w max ). 

While the specific choice of spin torque cutoff function 
V(p') (e.g. , we could have used a smooth cutoff, as op- 
posed to a sharp, Hcaviside cutoff) may be important 
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FIG. 4: Dissipative droplet frequency (dashed and solid 
curves) as a function of sustaining current from Eq. (|18[) . Pa- 
rameters are = 5.24 and v = 0. The dash-dotted line is 
the FMR frequency. 



in numerical applications^!, it has only a slight effect on 
the droplet sustaining current. For the Heaviside cutoff 
considered here, Eq. (fT5]l simplifies to 



r °°sin 2 e (p';u)p'dp' 



sm 2 Q (p';uj) 
o 1 + i/cos6 (/o'; oj) 



(20) 



p' dp 1 



Equation (|2T)|) reveals the explicit dependence of the sus- 
taining current on two key spin torque parameters: the 
contact radius p'„ > and the spin torque asymmetry 
< v < 1. We now investigate properties of the dissipa- 
tive droplet as a function of these two parameters. Fig- 
ures \5l a-d) show the dependence of cr™™, m™ ax , o; max , 
and the droplet radius p' dmp on the contact radius for var- 
ious choices of the spin torque asymmetry. The droplet 
radius is defined to be the radius at half maximum: 




1 



FIG. 5: Dissipative droplet properties for varying contact ra- 
dius p* and spin torque asymmetry v. (a) minimum sustain- 
ing current, (b) maximum frequency, (c) maximum m z at ori- 
gin, and (d) droplet radius (dotted line is p' dlop — p', , plotted 
for comparison). 



"Madron) = oC 1 



(21) 



Figure |5fa) shows that droplets excited by smaller con- 
tacts require larger sustaining currents. The dependence 
of a sus /a on p'^ and v for small contact radii can be made 
explicit by an asymptotic evaluation of the denominator 
in eq. ([2TJ|) giving 



2wfl 



v cos 



p?sin 2 e o (0;w) 



minimum droplet radius of about 2y/h^ — 1 l cx in dimen- 
sional units. Finally, spin torque asymmetry has only a 
small, perturbativc effect and does not substantially alter 
the droplet solution. 



+0(1), < pi « 1. 



(22) 

Therefore, a™^/a grows like l/pf for small contact radii 
in agreement with Fig. [5ja). Interestingly, we find that 
the droplet is constrained to have a frequency < uj < 
0.75 (Fig. EJb)), significantly below the FMR frequency 
of 1. Also, since m™ ax < (Fig. Etc)), the droplet al- 
ways has some region near its center that is partially 
inverted with respect to the surrounding magnetization. 
As shown in Fig. [S{d) , the droplet is well localized within 
the nanocontact; e.g. p' drop < p'*, when p'„ > 3. There is a 



2. Case v = 



We consider now the case without spin torque asym- 
metry (i.e. v = 0) in which case we can study the effect 
of varying the angle of the fixed layer Of on the droplet 
dynamics. In this section, we use the rescalings in Eq. ([5]) 
but do not use the transformation (|16[) . 

We solve Eqs. CS]) and (H^with v = 0, 6i(r',t') = 
6 s (r')sino;i', and $i(r',t') = $(r') + $ c (r') coswt', lead- 
ing to the following system of non-homogeneous equa- 
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tions: 

L<j,§ c +ojQ s = sin 0fWp') cos ©o 

a 

L e e s -Lu^ c = --sin6 l V(p'), 



(23) 
(24) 



— (uj + h' ) sin6 H cos8fV(p') sin© . 

a 

(25) 



Note that, in contrast to the case 9{ = 0, the overall 
phase <& no longer evolves linearly in time. Applying 
L® to Eq. (|23| and L$ to Eq. (|24|) gives the decoupled 
system: 

(LqL^, + lo 2 )$ c = — sm0f[LuV(p') 
a 

-L @ {V(-)cosQ (-)}(p')}, 

(L^Lq + uj 2 )Q s = -sm9t[L${V(-)}(p') 
a 

-(jV(p')cosG (p% 

These equations are always solvable if LqL§ + uj 2 and 
LqLq +uj 2 are strictly positive operators. One can show 
that L$ > 0, so that L^Lq +uj 2 > 0. One can also show 
that L e > — ?7q, where — t]q is the smallest eigenvalue 
of Lq and is strictly negative. Then, LqL$ > —T)q. 
We have verified by numerical computation of t]q that 
w 2 > t/q so that LqLq+uj 2 > as required, and Eqs. (|23|) 
and (|24[) are solvable. 

We are interested in the solvability condition for $ in 
Eq. ([25]), which is 



a sus (uj) = a(ui + h' Q ) secOf 



J °° sinQo(p';a;V dp' 
J^V(p') sin 2 Q {pf-u)(/dpf' 



(26) 

Note that the applied field appears only as a shift in the 
droplet frequency, as in the case 6{ = studied in the 
previous section. The expression (|26|) for the sustain- 
ing current agrees with the previously derived sustaining 
current in Eq. (|18[) when v = and #f = 0, as required. 
Thus, the dominant effect of rotating the fixed layer is to 
increase the sustaining current in proportion to sec 6[ . 



III. PHYSICAL PERTURBATIONS OF A 
DISSIPATIVE DROPLET 

So far, we have considered the dissipative droplet so- 
lution only for a simplified geometry where asymptotic 
methods can be applied. In these cases, the external field 
is both uniform and oriented perfectly perpendicular to 
the film plane. By making this geometrical simplification, 
we were able to factor out the contribution of the exter- 
nal field from the droplet solution. However, in a real 
point-contact system, the external field is neither uni- 
form nor perfectly perpendicular. In particular, the cur- 
rent flowing through the contact is an additional source 



of spatially inhomogeneous magnetic field, the Oersted 
field, and the applied uniform magnetic field in actual 
experiments is usually tilted away from the perpendic- 
ular axis. We employ micromagnetic simulations to in- 
vestigate how these physically important perturbations 
to the external field alter the ideal droplet solution. The 
numerical details used for our simulations are presented 
in Appendix [Bj 

We find that the combination of external field tilt and 
the Oersted field breaks the symmetry of the solution 
such that the droplet is no longer centered in the middle 
of the contact. As a result, the solution takes on a non- 
trivial inhomogeneous phase structure where the phase 
of the spin precession closer to the center of the contact 
precedes the phase further from the center of the con- 
tact. In addition, the spatial structure of the droplet is 
no longer perfectly circular. For some particular com- 
binations of simulation parameters, the droplet breaks 
away from the contact altogether and dissipates, a be- 
havior we call a drift instability. When this occurs, the 
droplet may maintain its form for many precession cy- 
cles, but it eventually decays, since outside the contact 
there is no spin torque excitation to balance damping. 



A. Oersted Field 

First, we consider the effect of the current induced 
Oersted field while keeping the applied field and fixed 
layer oriented almost normal to the film plane (canted 
by 0.00001° and 0.40°, resp.). The reason for this 
slightly asymmetric configuration is to test whether 
high-symmetry solutions are structurally stable to small 
changes in the system parameters. Such a configuration 
is experimentally possible, in principle. Our model for 
the Oersted field was presented in Ref. [l9| and takes the 
form 



-/oe(p)v, 



where 



/oe(p) = g oe (p) + 



2irM s r* 



p/p*, < p < p* 

p«/p, p* < p, 



(27) 



The function g{p) given in Ref. [jjl involves integrals of 
Bessel functions and depends on the geometry of the 
current-density distribution. The parameters defining g 
in Ref. |H are d = 1.67, z* = -0.925, and a = 2.92. 
The other, more dominant, term in Eq. (|27|) is the 
field generated by an infinitely long conducting wire. 
As an example, for the simulation in Fig. [51 we have 
max pe[0:Oo) \g oc (p)\ = 0.0081, while I/2nM s r* = 0.086, 
an order of magnitude difference. 

Figure [5] illustrates how the Oersted field changes the 
structure of the droplet. In contrast with the symmetric 
case in Fig. [3[ the azimuthal angle $ shows significant 
spatial variations. The droplet is also slightly shifted off- 
center. The strong phase variations are indicative of a 
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FIG. 6: Dissipative droplet in the presence of the Oersted field 
and with the applied field and fixed layer nearly perpendicular 
to the plane.. Parameter values are = 1.25, a = 0.03, 
p, = 5.24, h = 1.8, O = 0.00001°, 6t = 0.40°, v = 0.257, 
a = 0.196. 

tendency for the droplet to propagate^. In some cases, 
although not for the simulation in Fig. [6l the droplet 
breaks free from the nanocontact. This drift instability 
is discussed further in Sec. IIIIBI 

In Fig. [Jj we compare the numerically computed per- 
turbed droplet frequency as a function of current with 
the sustaining current from Eq. ([20)1 . We find that the 
droplet frequency is approximately shifted down by an 
overall amount of 0.012, but otherwise follows the same 
trend as the symmetric result. This behavior demon- 
strates that the analysis of the previous section, despite 
the necessary high symmetry restrictions, yields relevant, 
qualitative information about localized structures excited 
in a physically realizable nanocontact. 

As shown in Sec. IIV1 when the current exceeds a 
threshold value; e.g., the Slonczcwski critical current, 
a droplet can nucleate for sufficiently large anisotropy. 
The vertical line in Fig. [7J is the numerically computed 
threshold current. It differs from Slonczewski's result be- 
cause h oe ^ 0, 6> > 0, and 6{ > 0. This threshold for 
droplet nucleation suggests a hysteretic effect that will 
be discussed in Sec. IIV CI 

Our micromagnetic simulations also show that the per- 
turbed droplet, for certain parameter choices (e.g., the 
crosses in Fig. [7]), undergoes a drift instability. This be- 
havior will now be investigated further. 

B. Canted Applied Field, Fixed Layer, and 
Oersted Field 

In this section, we investigate the combined effects of 
the Oersted field as well as canting of the applied field 
and fixed layer. Figure [8] is a time sequence showing 
the evolution of a strongly perturbed droplet over one 
precessional period. In contrast to the nearly symmetric 
configuration of Fig. [51 where the droplet was slightly 
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FIG. 7: Dots: perturbed droplet frequency as a function of 
current with the Oersted field computed from micromagnetic 
simulations. Crosses: perturbed droplets that undergo a drift 
instability; the frequency is calculated before the instabil- 
ity manifests. Solid curve: droplet frequency as a function 
of the sustaining current from Eq. (|20[) . Dash-dotted: the 
Zeeman (ho) and FMR (ho + h^ — 1) frequencies. Triangle: 
the Slonczewski critical current and onset frequency for high 
symmetry^ (see Appendix [A|l . Solid vertical line: numerically 
computed threshold current in the presence of the Oersted 
field and the canted fixed layer. Parameter values are the 
same as in Fig. [6] 



shifted to the left, the droplet is slightly shifted down, 
toward the region of lower in-plane field. 

As the system parameters are changed, the shifting of 
the droplet center can be large enough to actually dis- 
lodge the droplet from the nanocontact. An example of 
this drift instability is shown in the panels of Fig. [SJ The 
current was taken to be less than the current for Fig. [8j 
which did not experience a drift instability. A number of 
precessional periods pass before the droplet breaks free. 
Once free, it propagates, but because it no longer sat- 
isfies the required balancing condition between damping 
and spin torque, it loses amplitude and decays. Once the 
droplet has drifted outside of the nanocontact, a new one 
is formed if the nucleation conditions arc satisfied (sec 
Sec. IIV[) . Wc have also observed droplets that drift sev- 
eral nanocontact diameters before decaying, i.e. the cen- 
tral magnetization lifts up so that min rg K 2 w«(r, t) > 0. 
The manifold in parameter space in which the droplet 
manifests a drift instability appears to be complicated. 
Nevertheless, wc readily find parameter regimes where 
the droplet apparently does not experience the drift in- 
stability, as in Figs. [6j [Jj and[8j 

Notice that the droplet propagates down, in the — y di- 
rection. Recall that the canting direction of the applied 
field is along x. Due to the symmetry of the Oersted 
field, the direction of propagation of the drifting droplet 
appears to track the azimuthal angle of the applied field 
minus 90°, if 9q is sufficiently large. For example, if the 
applied field is canted along y, then the droplet will drift 
along x if unstable. This can be understood as a conse- 
quence of magnetostatic interactions between the effec- 
tive dipolc moment of the droplet and the field gradi- 
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FIG. 8: Time sequence of a strongly perturbed droplet over one period of precession in the presence of a canted applied field, 
canted fixed layer, and Oersted field. Parameter values are = 1.25, a = 0.03, p* = 5.24, ho = 1.8, 9o = 5°, 8f — 31.4°, 
v = 0.257, a — 0.189. The time of the initial panel here and in Fig. [9] is set to for comparison. The simulation actually began 
earlier. 



ent associated with the Oersted field. Given the canting 
of the applied field, the effective dipole moment of the 
droplet acquires an in-plane component that is drawn to 
the edge of the contact where the Oersted field is also 
in the — x direction such that the Oersted field gradient 
acts to trap the droplet. In the case of Fig. [51 where the 
applied field is canted only 0.00001°, the droplet is ob- 
served to drift to the left rather than down. Therefore, 
the strength of the Oersted field and the in-plane com- 
ponent of the applied field have a strong effect on the 
existence and dynamics of a drift instability. 



IV. NUCLEATION OF A DISSIPATIVE 
DROPLET 

Figure [10] shows the birth of the droplet pictured in 
Fig. [5] starting from a state pointing uniformly in the z 
direction. For sufficiently large perpendicular anisotropy 
/ik and current, the magnetization inside the nanocon- 
tact reverses and nucleates a droplet. In this section we 
show that the reversal mechanism is caused by an insta- 
bility of small amplitude waves. We will study Eq. (TTJ) in 
the weakly nonlinear regime and find that the small am- 
plitude Slonczewski mode that exists near the threshold 
for the onset of dynamics is stable/unstable depending 
on whether is less than/greater than a critical value 
h". We find that > 1 due to exchange effects and 
converges to 1 as the contact size is increased. The value 
of h" is important for the possible experimental obser- 
vation of a dissipative droplet. 



A. Stability Analysis of a Macrospin 

Before studying the PDE ((T|), we consider the 
macrospin model where spatial variation is neglected. 
The stability analysis for this model is suggestive and 
mathematically simpler. However, we find that the crit- 
ical anisotropy field in the macrospin case /i£ r ' m satis- 
fies /i"' m < 1 and depends strongly on the applied field, 
which differs from the result obtained by analyzing the 
full PDE model, where we find hf* > 1 with weak applied 



field dependence. 

We consider Eqs. and (0| neglecting all spatial vari- 
ation and inter-layer dipole coupling in the symmetric 
regime 8q = and 6f = 0; i.e., we are assuming that the 
initial condition for the system is in the parallel state. In 
this case, the equation for O is decoupled from $, so wc 
can just study the scalar, first-order ODE 



6 = — asinO[cosO(/ik — 1) + h ] 



■ sin O 



1 + v cos 6 



(28) 



0, we 



By linearizing Eq. (|28j) about the equilibrium O 
find that it becomes unstable when 

a > cr = a(l + v){h Q + h k - 1), 

in agreement with previous numerical and mathematical 
analyses of this system contained in Ref. H^. We seek 
a periodic equilibrium solution O(t) = O c just above 
threshold by taking 



Then, wc have 



Oe 



(T + £, < £ <C (T . 



2ae 



1/2 



.l-/i k + i/[2(l-/i k )-/i ]. 
This solution exists (is real valued) as long as 



ftk < IK = 1 - 



1 + 2^ 



(29) 



Furthermore, one can show that this equilibrium is sta- 
ble. Therefore, when < h"' m , the equilibrium O = 
undergoes a supercritical Hopf bifurcation as the current 
a is increased beyond (To- When > /i k r ' m , there is no 
periodic solution, and the system switches from O = to 
the fully reversed state o = ir when a exceeds ctq. 



B. Stability Analysis of the Micromagnetic System 

We consider Eq. (fT]) in the symmetric regime with 8q = 
0, h oc = 0, 8{ = 0, and the substitutions 



u = m x + im v , m z 



10 



t = o 



t = 195 



t = 205 



t = 214 



t = 223 



t = 232 




n 1 



FIG. 9: Time sequence showing the droplet drift instability for the same parameter values as in Fig.[8]but with smaller current 
a — 0.121. To facilitate visualization, the length of the in-plane magnetization vectors is normalized to the largest value in 
each frame. 
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FIG. 10: Birth of a dissipative droplet soliton for the current a — 0.186, above the Slonczewski critical current er s = 0.160. 
Parameter values are the same as those in Fig. [8] 



where | u | <C 1. Then u approximately satisfies a complex 
Nonlinear Schrodinger type equation 



du 
l dt 



(1 + ia)\7 2 u — (ho + — l)u — ia(h a + h k — l)u 



(uV 2 |u| 2 - |u| 2 V 2 u) + ia\Wu\ 2 u 
aV(p) 



a(h + 2h k - 2) - 



{1 + uf 



u\ u. 



(30) 



Similar to the macrospin case, when a is increased past 
a threshold value, the zero solution becomes unstable. 
The threshold, critical current and onset frequency were 
found by Slonczewski as a solution to a linear eigenvalue 
problem^ (see Eq. (|A1I) ) . Incorporating weak nonlin- 
ear effects, Ref. [24| showed that a small amplitude pe- 
riodic solution exists as a modulation of the Slonczewski 
mode for h k = and v = 0. This time-periodic, weakly 
nonlinear mode was found to be numerically stable. In 
this section and Appendix [XI we extend these results to 
< /ik < and < v < 1 , where is defined through 
Im(£(/i")) = and £ is given by Eq. (|A4|) . This behavior 
is analogous to the supercritical Hopf bifurcation for the 
macrospin, studied in the previous section. 

In this section, we show that weakly nonlinear mod- 
ulations of the Slonczewski mode are no longer stable 
when hk > > 1. This behavior is analogous to the 
switching exhibited by the macrospin for hy_ > h^' m and 
a above threshold. While we are unable to analytically 



follow the dynamics of the instability, numerical simula- 
tions such as the droplet birth sequence shown in Fig. [10] 
demonstrate that the magnetization reverses inside the 
nanocontact and develops into a dissipative droplet. 

For the stability analysis, we seek a multiple scales 
solution representing a modulation of the Slonczewski 
mode in the form 



u(r,t) 



e i{ - ho+hk - l)t e iUst [eA(T)f(p) 



+ e 3 Ul (r,T) + 



(31) 



£ 2 (7i 



where < e <C 1 is the amplitude of the mode at the 
origin, which is modulated by A(T) with T — e 2 t the 
"slow" time, w s is the frequency of the Slonczewski mode 
f(p) with threshold current er s , and <r\ represents a de- 
viation from a s . The explicit form for / and the implicit 
equations for ui s and <r s are given in Appendix |A"1 

Invoking a solvability condition at 0(e 3 ), Eq. (|A3[) 
gives the nonlinear amplitude equation 



dA 



0"! 



1 



■ (A + tlAfA, 



(32) 



with complex linear and nonlinear coefficients £ and £. 
There is a time-periodic solution of Eq. (|3^|) for a specific 
choice of o\ € M 



-(1+^) 



A(T) 
Re(C) ' 



MO 



MO M0 



Re(C) 



(33) 
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FIG. 11: Modulation parameter Im(£) as a function of h^. 
When Im(£) > 0, the weakly nonlinear Slonczewski mode 
is modulationally unstable. Parameter values are a = 0.03, 
p, = 5.24, h = 1.8, v = 0.26. 
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FIG. 12: The critical anisotropy field ftjf as a function of con- 
tact radius p, for various spin torque asymmetries v. Other 
parameters are a = 0.03 and ho — 1.8. 



is small. This is confirmed by numerical calculations of 
the h" dependence on ho > for the values of p* and v 
plotted in Fig. Q21 We find that h" varies by at most 3 % 
for ho € [0,1.8] with j^h"(p*,v,ho) < 0. Importantly, 
the lower bound for /i£ r is preserved: h™(p*,v,ho) > 1 
when ho > 0. This behavior stands in stark contrast 
to the macrospin result for h"' m < 1 in Eq. (f2T)j) that 
strongly depends on the applied field. 

From Eq. p3|) . we see that o\ is negative when 
Im(£)/Re(C) > 0. Numerically, we find that Re(C) > 0, 
so that a i < when h^ > h", which corresponds pre- 
cisely to the instability criterion. Furthermore, when 
cr i < 0, any finite-amplitude excitation will tend to grow, 
thus the time-periodic solution is unstable. This behav- 
ior is similar to the macrospin model discussed in the 
previous section where, for hk > h", a small amplitude 
periodic solution did not exist. 

Numerical simulations of Eq. ([T]) confirm the foregoing 
analysis, even in the non-symmetric cases with 9o > 0, 
6f > 0, and h oc ^ 0. There is a critical value of the cur- 
rent above which large amplitude dynamics ensue. When 
/ik is above h" , we find that the magnetization inside the 
nanocontact reverses to form a localized, coherently pro- 
cessing, fully nonlinear magnetic solitary wave. We iden- 
tify this solitary wave with the dissipativc droplet soliton 
found in the asymptotic analysis of Sec. [TTJ 

Due to the numerically robust formation of the droplet 
for a variety of initial data and across a large parameter 
regime, we view it as a global attractor. As long as the 
current is above threshold, a droplet is observed to nu- 
cleate. 



This solution represents the nonlinear frequency shift 
£ 2 w s ,i to the Slonczewski frequency uj s . The stability 
analysis in appendix [3] shows that the solution (|3"3"]l is 
unstable when Im(£) > 0. An explicit formula for £ is 
given in Eq. (|A4[) . We evaluate the integrals numerically 
and plot Im(£) as a function of hy for specific param- 
eter values in Fig. [TT] There is a critical value of the 
anisotropy field h" satisfying 

Mtfftk )] = 0, (34) 

above which the weakly nonlinear Slonczewski mode is 
modulationally unstable. In other words, weak modula- 
tions of the Slonczewski mode will grow exponentially in 
time when hu > h"' . Figure [T2"l shows the dependence 
of h£ = h^{p*,v,ho), computed numerically by solving 
Eq. (|3"4"|) with Eq. (|A4[) . as the contact radius and spin 
torque asymmetry are varied and ho = 1.8. We see that 
h" is strictly greater than one and that larger perpen- 
dicular anisotropy is required to enable the nucleation 
of a droplet in smaller nanocontact devices. For a given 
contact size, larger spin torque asymmetry permits nu- 
cleation of a droplet at smaller values of h^. 

We have also investigated the dependence of h" on the 
applied field magnitude ho- Since ho appears in £ only 
multiplied by a, we have gf-^kX/ 9 *' v i ^o) = 0(a), which 



C. Hysteresis 

We have shown in this section that a dissipativc droplet 
will form when the current exceeds the threshold for in- 
stability of small amplitude waves inside the contact. In 
the high symmetry case of Eq. @ and Of = 0, the thresh- 
old is the Slonczewski critical current a s , which is plotted 
in Fig. [7] (triangle) along with the numerically calculated 
threshold current from micromagnctic simulations that 
incorporate the Oersted field and a canted fixed layer 
(cr = 0.156, vertical line). We observe that the mini- 
mum droplet sustaining current is below the threshold 
for droplet nucleation. This suggests a hysteretic effect, 
whereby a droplet can be nucleated at a current above 
threshold and remains stable when followed by an adia- 
batic decrease of the current below threshold. We have 
performed this experiment numerically for the perturbed 
droplets of Fig. [7] We use an already nucleated droplet 
(at a = 0.211) as the initial condition for a new simu- 
lation with cr = 0.151, below threshold. We find that 
sufficiently close to, but below, threshold, the droplet's 
frequency slightly increases but remains stable. For fur- 
ther decrease of the current to a = 0.136, however, the 
droplet undergoes a drift instability similar to behavior 
shown in Fig. [51 Since the current is below threshold a 
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new droplet does not form. 



V. DISCUSSION 

The Slonczewski form for the spin torque term we are 
considering here, Eq. ([T]), was derived strictly for the case 
of spatially uniform magnetization distributions^. Non- 
local generalizations of the spin torque term to nonuni- 
form magnetization distributions have been derived in 
the small amplitude regime^ and applied to single layer 
nanocontact simulations^. Because the spin torque and 
damping terms are treated as perturbations giving rise to 
a dissipative droplet solution, we expect that the specific 
form for the spin torque term, and the damping term 
for that matter, will not yield qualitative changes in the 
structure of the droplet. The asymptotic analysis pre- 
sented here is applicable to an arbitrary spin torque term 
as long as its net effect is to oppose the inherent damping 
in the system. 

Due to the symmetries in the problem when 9 = Of s» 
0, one might expect the Oersted field to prefer the exci- 
tation of a localized vortex structure; e.g., a topological 
soliton^. We observe no such excitations in our simula- 
tions. Two-dimensional, conservative topological solitons 
are stable but they have higher energy than the conserva- 
tive non-topological soliton^, providing one possible ex- 
planation for this behavior. Another possible explanation 
lies in the form of the spin torque term. In contrast to the 
local form for the spin torque considered here, single-layer 
nanocontact simulations that incorporate a nonlocal spin 
torque have demonstrated the formation and numerical 
stability of a precessional vortex in high-symmetry con- 
figurations without perpendicular anisotropyi^. A more 
realistic model of the spin torque may favor a topological 
soliton in certain cases. 

Unfortunately, we as yet have no detailed physical un- 
derstanding of either the dissipative droplet asymmetry 
or the drift instability in the presence of the combined 
effects of an Oersted field and a tilted applied field. In 
part, our lack of understanding stems from the limited 
theory for the spatial propagation of the drople t 22 ' 26 . Ex- 
tension of the theory presented here to cases where the 
droplet is accelerated by field gradients and other such 
forces may provide more physical insight to aid in un- 
derstanding why the droplet is subjected to displacing 
forces. 

The physical appearance of the droplet is reminiscent 
of the magnetic bubble that was once the subject of in- 
tense investigation as a possible alternative to fcrrite- 
core computer memoryi 6 -. Indeed, we can identify the 
zero- frequency droplet as a topologically trivial magnetic 
bubble^ with a winding number of zero, though in the 
theory presented here, the droplet frequency approaches 
zero only in the limit of infinitely large diameter. How- 
ever, we expect that the droplet will converge to the bub- 
ble structure at a finite diameter, contingent on the inclu- 
sion of non-local magnetostatic energy in the calculation; 



even when considering free layer films as thin as 3 nm, 
such a term will eventually overcome the positive con- 
tribution of exchange energy to the oscillator frequency, 
stabilizing the droplet as a static structure. In this sense, 
then, we can think of the conservative droplet soliton, 
a delicate balance between exchange and perpendicular 
anisotropy, as a dynamically collapsing bubble, and the 
dissipative droplet soliton as an imminently collapsing 
bubble that is critically stabilized by the localized injec- 
tion of spin torque. 

VI. CONCLUSION 

We have derived equations for a dissipative droplet 
soliton through an asymptotic analysis of the Landau- 
Lifshitz-Slonczewski equation for a point-contact spin 
torque oscillator with perpendicular anisotropy in the 
free layer. The droplet soliton is a localized, dynamic, 
solitary wave solution consisting of partially reversed 
magnetization directly under the contact and a zone of 
large amplitude precession in a region bounding the re- 
versed magnetization. The diameter of the processing 
boundary is approximately the contact diameter. The 
droplet frequency is always strictly less than the ferro- 
magnetic resonance frequency for the film, and is also 
a monotonically decreasing function of droplet diameter. 
The balance between spin torque and damping required 
to sustain the droplet determines the relation between 
bias current and frequency. The instability of small am- 
plitude solutions that leads to formation of the reversed 
domain in the droplet requires a minimum perpendicular 
anisotropy that is a function of the contact radius and 
spin torque asymmetry. 

Some of the unique, identifying properties of the dis- 
sipative droplet that could be observed experimentally 
include 

• The droplet 's frequency is well below the ferromag- 
netic resonance frequency. 

• Sufficiently far from the minimum sustaining cur- 
rent, the droplet's frequency has a weak depen- 
dence on current. 

• The droplet may manifest a drift instability, which 
would reveal itself as a transitory cessation in ac 
oscillations until the droplet again forms. 

• Hysteresis in current is expected, unless a drift in- 
stability results in a finite droplet lifetime. 

We have investigated the nucleation and stability of 
the dissipative droplet soliton through numerical simu- 
lations. We find that droplet formation begins once the 
current in the point-contact is sufficient to instigate the 
small amplitude Slonczewski mode, characterized by spin 
waves that radiate away from the point-contact. For suf- 
ficiently strong perpendicular anisotropy, this mode is 
subject to a modulational instability and rapidly evolves 
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into the reversed magnetization profile of the droplet soli- 
ton. We find that the droplet is stable in certain param- 
eter regimes with regard to the inhomogeneous Oersted 
field and to variations in spin torque asymmetry and ap- 
plied field angle. Finally, the droplet is subject to a drift 
instability that is a complicated function of the parame- 
ters employed in this theory. 



has been chosen so that Im(fc ) > 0, and / experiences 
exponential decay due to material damping a > 0. The 
decay length is weak compared to the contact radius. The 
two real eigenvalues lj s and a s are determined by solving 
the complex valued transcendental equation 

hH^ikoP^Mhp.) = k o H[ l \k oP *)J {h P *), 



Appendix A: Modulational Instability of 
Slonczewski Mode 



which results from continuity of the first derivative of /. 

Continuing the asymptotic analysis to the next order 
gives the following nonhomogeneous equation 



Here we provide the details of our stability analysis for 
small amplitude, modulated waves excited in a nanocon- 
tact. 

Inserting the ansatz (|5T|) into Eq. (|3T)|) and consider- 
ing the leading order behavior in e gives the linear Slon- 
czewski eigenmode f(p) satisfying^ 



L f = (1 + ia)(f" + -/') - ia(h + K- 1)/ 
P 

1 + V 



(Al) 



Slonczewski considered both u s and a s as eigenvalues for 
this equation and found the C 1 [0, oo) solution 

,, n. / M k iP), < P < P* = Jojhp*) 
nP) ~\cH^\k oP \ p*<p ' ° ffW^p.)' 

where Jo is a Bessel function and Hq is a Hankel func- 
tion. The inner and outer wavenumbers fci, k are 



h = 



k Q — 



uj s - i[a(h + hk - 1) - er s /(l + v)] 



1 + ia 



uj s - ia(h + h k - 1) 



1 + ia 




v s V(p) 



(1 + uf 



l/l 2 / 



|/| 2 [/"+ V] + 2ia\f\ 2 f 



(A2) 



Since the kernel of the adjoint of Lo is spanned by /* (* 
denotes complex conjugation), we invoke the solvability 
condition of Eq. (|A"2j) 



f(p)R 1 (p)pdp = 0, 



(A3) 



to determine the dynamical equation for A(T) in 
Eq. ([52)1 . The complex valued linear and nonlinear coef- 
ficients C and £ are 



[Jo(fcijO*)M(fciP*)]' 
1 - [fci/fco] 2 



Since \f(p) \ ~ Ce" Im(feo)p / v /p for p > p*, the sign of k Q and 



(i 



(1 + 1/)= 



[h k - 1 + ia(h + h k - 1)] / \f\ 2 fpdp - i - / \f\ 2 .fpdp - (1 - 2ia) I \f'\ 2 fpdp 



(A4) 



plJi{kp*y\i- (fe/fco) 2 ] 

i 

We now perform a linear stability analysis of the time- which experiences exponential growth when Im(£) > 0. 
periodic solution A(T) = e lUs > lT in Eq. (|33| by inserting 

A(T) = e i£ " s ' lT (l + v + iw), v,w&R, \v\,\w\<l, 

Appendix B: Micromagnetic Computational 
into Eq. (|32|1 . Keeping only the terms linear in v and w Modeling 
gives the simple, decoupled dynamical system 

dv dw The numerical method we have used to solve Eq. ([T]) 

—j^ = 2Im(£)f, -= = — 2Re(£)u, j g s i m il ar to the method tersely presented in Ref. [l9l but 
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adapted to a trilayer nanocontact. In this Appendix, we 
describe our method in detail. In brief, we use a polar 
coordinate grid and cylindrical magnetization basis. An- 
gular derivatives are computed by use of a pseudospec- 
tral, Fourier method while radial derivatives on a nonuni- 
form grid are computed by use of finite differences. For- 
tune-stepping, we use an explicit, hybrid 2 nd /3 rd -order 
Runge-Kutta time stepper with error control, ode23 from 
Matlab®. 

The polar coordinate system is a particularly efficient 
and accurate choice for nanocontact simulations. The 
discretization we use is non- uniform in radius ("inner" 
and "outer" grids) 



tanh | - — - I 



Pn = q(n) = p* j | (dpout - dpin) 
+ dp in X d£, 



(Bl) 



where n = 1, . . . , N, and n, w are parameters determin- 
ing the location and width of the smooth change from 
the fine inner grid spacing dp ln to the coarser outer grid 
spacing dp out > dp ln . We typically have p\ ps dp ln . For 
the angular, pseudospectral discretization, we expand the 
cylindrical magnetization basis in a truncated Fourier se- 
ries with M Fourier modes 



m(p n ,<p,t) 

M/2-1 



E 

fc=-Af/2 



Jk<p 



» 



(B2) 



The transformation from an cquispaced angular grid 
tpk = —n + (k — l)2n/M with magnetization coefficients 



L n,k 



(t) evaluated at the discrete grid (p n ,'Pk), to the 
Fourier representation in Eq. (|F32|) is achieved by use of 
the fast Fourier transform. 

The advantage of this discretization is that we can 
solve on a uniform computational grid while the physical 
grid is clustered in and around the nanocontact, where 
the majority of the fine scale dynamics occur. The outer 
grid supports the propagation of spin waves of the appro- 
priate wavelength away from the nanocontact. In order 
to simulate an infinite domain, we choose a finite domain 
large enough, L = piq ~ 30p*, so that spurious waves are 
naturally damped to a negligible amplitude. Then their 
reflection off the boundary does not affect the strongly 
localized dynamics near the nanocontact. 

Computing the vector Laplacian V 2 m is the compu- 
tationally intensive portion of the algorithm. Angular 
derivatives are approximated by multiplication in Fourier 



space 



dip 



ik. The approximation of radial derivatives 



using finite differences requires some care, especially near 
the origin. We discuss the details now. 

Due to the non- uniform radial grid (|B1[) , radial deriva- 
tives in computational space require appropriate factors 



of q and its derivatives. For example 



d_ 

dp 



\_d_ 

q' dn 



dp 2 



r 



dn 2 



<1 



q 



dn 



The radial derivatives in computational space are approx- 
imated using 6 order finite differences. The boundary 
condition at the artificial boundary p = L for Eq. ([T]) 

is the Neumann condition — — (t) = 0, which is im- 
plemented by use of a standard ghost point method and 
one-sided differences. For radial derivatives near the ori- 
gin at n = 1, 2, we define symmetry conditions to obtain 



;(P>¥>>- 



(t) by taking ip — > ip + ir: 



(B3) 



where we have used p — > —p and ip —> —ip when ip —> 
p+n. Therefore, fn^'^\t) are even/odd functions of n as 

k is odd/even while fh^\(t) is even/odd as k is even/odd. 
At p = 0, we take (see Ref. [H) 



"0,k 



(t)=0, \k\^l, 



l 0,k 



(t) =0, k £ 0, 



dm 



0,±1 



(B4) 



dp 



(t) = 0, 



dp 



= 0. 



The derivative conditions in Eq. (|B4j) are approximated 
by use of one-sided differences to obtain an estimate of 
rfbQ jfj (t) and mQg(t). The symmetries in Eq. (|B3[) and 



the conditions in Eq. (|B4|) enable the use of centered fi- 
nite differences, even near the origin. An explicit, Runge- 
Kutta 2 nd /3 rd -order time-stepping method is used to ad- 
vance the discrctized version of equation ([!} forward in 
time while renormalizing the magnitude of m after every 
time-step to preserve the constraint |m| = 1. To avoid 
severe time-step restrictions due to the small grid spacing 
near the origin (w 2np n /M), we apply a smooth, radial 
grid dependent angular mask that effectively reduces the 
number of angular modes at p n from M to M n = 2ir/k n : 



Qn(k) = — I — tanh [ — — L 
y K ' 2 2 \ Ak 



The mask's parameters are the radial grid dependent 
wavenumber cutoffs k n , n = 1,...,N and the cutoff 
width A A;. Application of the mask at every time step 
filters out numerically induced small wavelengths near 
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the origin^. The cutoffs k n are chosen so that the ap- 
proximate grid spacing 2ir p n /M n ~ dp- ln , hence the grid 
near the origin has an effective spacing of dp- ln . The mask 
applied to rh n k (t) takes the form 

gW{fnW(t)} = g n {k)rh { :\(t). 

I. J n,k 

Care must be taken when applying the mask to the in- 
plane Fourier coefficients rh^'^ (t) because p and (p de- 
pend on the grid location. We use 

{m^(t)} n k = \{[g n {k 1) + g n (k + l)]m$ (i) 

+ i[-g n (k-l)+g n (k + l))m%l(t)}, 
{m<*> (*)} n k = \{i[9n{k 1) - 5 „(fc + (i) 

+ [9n(k - l)+g n (k + l)]fh^l(t)\. 

* Electronic address: mahoefer@ncsu.edu 
t Contribution of the U.S. Government, not subject to copy- 
right. 
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